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Abstract. Skew density matrices can be diagonalized to yield probability interpre- 
tation. The power-counting prediction of perturbative QCD is found consistent with 
recent CEBAF data on F 2 {Q 2 )/F 1 (Q 2 ). 
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Perturbative QCD can be applied to hard exclusive processes in many ways. Use 



of "skew" or "off-diagonal" parton distributions [f] has generated attention. Skew 
distributions are matrix elements sharing features with density matrices. The diag- 
onal elements of a positive density matrix are interpreted as probabilities. Density 



matrices depend on the basis: probability in one frame will appear to be an inter- 
ference of amplitudes in another frame. This occurs with quark spin amplitudes 



| [2], where transverse polarization ( "transversity" ) distributions appear as interfer- 

ence in the helicity basis, and helicity distributions appear as interference in the 
transverse polarization basis. 

Here we find a probabilistic interpretation for skew distributions involved in the 
proton electromagnetic form factors. Then a simple, kinematic angular momentum 
argument resolve a recent puzzle from data for the electromagnetic form factors 
Fi(Q 2 ),F 2 (Q 2 ). For about 30 years it was thought that F 2 (Q 2 ) / F^Q 2 ) ~ 1/Q 2 
was a prediction of hard scattering, and in particular, of pQCD. We made the same 
prediction on the basis of skew distributions in 1993 [3]. Recent CEBAF data [4] 

shows that Fi/F\ ~ 1/^J\Q 2 \. Our prediction missed some instructive points, and 
the CEBAF result is expected, if proper kinematics is simply taken into account. 

Diagonalization: Consider 2 — > 2 quark-proton scattering with momentum trans- 
fer A M , A 2 < 0. The proton is described by states \p — A/2, s > coming "in" and 
\p + A/2, s'> going "out". The scattering matrix element is 

$^(p; k, k', A)8(k +p — k' — p — A) 

= Jd A zd A z' exp{-ikz + ik' z') <p - A/2, s\ifp(z)^ a (z')\p + A/2, s'> . 



Here a,/3 are Dirac indices of the quark fields ip, and the in- and out-quarks have 
momenta k, k' with k = x(p — A/2) + k T , k' = x'(p + A/2) + k' T . A time-ordering 
symbol is dropped, as only one time ordering contributes: see Diehl and Gousset [1] 
Now $(A) can be decomposed into terms symmetric and antisymmetric in A; the 
symmetric part is Hermitian. This leads to a density matrix we can diagonalize. 

To diagonalize we make a series of coordinate transformations. Make a Lorentz 
transformation to the frame: p = (p + , 0, m 2 /2p + ); A = (0, A T , — t/2p + ); A T = 
—t. Now A M becomes entirely transverse as p + — > oo. The partons have 4-vectors 
k^ = xp^ + k T , x(p + A) u + k% , yielding 

x = x'\ k' T = kr + At (I — x). 

The matrix element is now diagonal in x. This might seem impossible, because 
the x,x' dependence of skew distributions is thought to be invariant: but x is 
not Lorentz- invariant sideways. The convolution in kT is diagonalized by conjugate 
transverse spatial coordinates b. Then jd 2 b s'; b/(l — x),x)e lb ' A can be inverted 
by Fourier transform to find the integrand, diagonal in everything but spin, which 
also can be made diagonal by familiar helicity or transverse bases. 

The choice of frames and diagonalization was used in an independent early in- 
troduction of off-diagonal distributions [3]. Due to a kinematics goof we missed the 
factor 1/(1 — x). A review by Brodsky and Lepage [5] was useful. 

Interpretation: The electromagnetic form factors (with — Q 2 = A T ) are found by 
following the Feynman rules, which includes multiplication by the quark charge e q , 
tracing $ with 7^, and doing the integral / dx d 2 b. At large — Q 2 this is dominated 
by b 2 ~ 1/|Q 2 | by Fourier analysis. The "short distance" implied by large Q 
is far more general than the more problematic "quark-counting" argument. We 
simply have one quark located by the hard momentum Q, while the asymptotic 
short distance theory [5] assumes that all possible Fock states are separated by 
asymptotically short distance. We remain within the framework of pQCD, of 
course, while choosing a more general factorization method. 

Since the form factor Fi(Q 2 ) is known from data, we invert the Fourier transform 
to solve for a positive-definite diagonal element of the density matrix, namely a 
probability: 

P(A, A; b) = \ [ d 2 Q T e'^F l {Q 2 ). (1) 
e q J 

Positivity is assured because F 1 (Q 2 ) is monotonically decreasing. The quantity 
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P(X, A; b) is on a similar footing to the usual parton distributions. Indeed, the usual 
parton distributions are functions of x integrated over fcy, while P(b) depends on 
b = \b\ and has been integrated over x. 

A fit to the Fourier transform of Fi(Q 2 ) has been performed. A profile of the 
transverse (1/(1 — x)-weighted) probability to find a quark (mostly of the up-type) 
was shown at the meeting. We expect that this invariant quantity will be useful in 
many studies involving the transverse coordinate. 



Orbital Angular Momentum: Quark orbital angular momentum is a fascinating 
subject of great interest in the proton spin puzzle. In the case at hand, we are lucky 
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to have two quarks evaluated at the same x, 6 points, so that difficulties of gauge 
invariance are minimal. Indeed, P(b) is gauge invariant by definition in terms of 
observable quantities. 

We expand the operators in terms of solutions to V 2 i/> = 0. By usual methods 
[2] the operators are evaluated inside the proton state, letting the correlation $ be 
expressed by c-number "wave functions". Partons are below threshold in a form 
factor, so we need an expansion for orbital angular momentum for spacelike k: 

ij{z-,b,z + = 0) = J2[dx i> n {x)I n {b)e- in «e ixp+z - . (2) 

Here e~ m< i> are light-cone 5*0(2) orbital angular momentum basis functions; I n (b) 
are modified Bessel functions more usually seen as J n {b) for a timelike basis. 1 We 
are interested in the F 2 form factor, associated with ia )lv Q v /2m, which represents 
proton helicity-flip at large Q 2 . Since we cannot change the helicity of a quark with 
a hard scattering (the near-perfect chiral symmetry of pQCD), the proton can only 
flip its spin to make F 2 by transferring a unit of orbital angular momentum in the 
quark [3]. (Hoodbhoy and Ji [6] subsequently verified the same result.) 
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Note we are not attempting here to derive the functional dependence of P(X, A; b): 
in our approach this comes from data. Our approach is to relate each power of b 
or angular momentum to further suppression by powers of 1/Q. The argument 
for short distance is kinematically compelling here (if controversial in the quark- 
counting method), so for large Q the Fourier transform is dominated by I (b) ~ b° 
if this channel is allowed. But I (b) represents the s-wave component, with zero 
angular momentum, and so this channel is open only to the helicity non-flip, namely 

Ft ~ 5 XX , J dx d 2 b Mx)TM%) [h (6/(1 " x))} 2 exp(iQ T • 6). (3) 

Here T represents the necessary Dirac matrices. 

Conversely, the only possibility for the helicity-flip F 2 is to use powers of b to 
conserve the angular momentum, and suffer the corresponding power-suppression 
in Q. On the basis of this power counting, it was reasoned [3] that by angular 
momentum selection rules, the integrals over b would vanish unless two representa- 
tions of the same angular momentum matched up, giving the previous prediction 
F 2 {Q 2 ) I F ± (Q 2 ) ~ 1/Q 2 for Q 2 > GeV 2 . Now consulting the formulas this is simply 
not true. There is a factor of e l ® T ' bT carrying angular momentum: 

F 2 ~ 6 X -v J dx d 2 b $ (a;)Ityi(a;)/o (6/(1 - x)) h (6/(1 - x)) j\QT\*™{*-<l>Q) + cc . 

(4) 

^ If one is concerned about b — > oo, then Green functions can be expanded in scries of 
7„(&<)i ; i'„(6>) where &<(&>) is the smaller (larger) of two b arguments. We only need short 
distance. 



Physically, the probe Q breaks the rotational symmetry of the problem. 

We reiterate that this analysis is entirely within the context of pQCD. In pQCD 
one takes some matrix elements from the data, and makes predictions for others. 
What can we predict here? From the power-counting cited, we have 

f 2 (Q 2 ) = <Mx)Mx)> (5) 

Fi{Q 2 ) y/\Q^\<bip (x)M^> 

where the braces represent the integrals. The fact that (Q/GeV)F 2 (Q 2 )/Fi(Q 2 ) is 
not far from unity in the CEBAF data indicates that the proton wave functions for 
quark angular momenta 1 and are not too different in magnitude. Constituent 
quark, or non-relativistic models are ruled out, but those models were never ca- 
pable of capturing the Fock space description of the skew distribution. It will be 
interesting to continue these studies in the context of the larger proton spin puzzle. 
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